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We analyse recent results on charged particle pseudo-rapidity densities from RHIC 
in the framework of the Dual String Model, in particular when including string fu- 
sion. The model, in a simple way, agrees with all the existing data and is consistent 
with the presence of the percolation transition to the Quark-Gluon Plasma already 
at the CERN-SPS. It leads to strict saturation of the particle (pseudo-)rapidity den- 
sity, normalised to the number of participant nucleons, as that number increases. 
Asymptotically, as V* — * 00 1 with the number of participants fixed, this density 
approaches again nucleon-nucleon density. A comparison with recent WA98 data 
is presented. 



The dependence of measurable quantities like charged particle density, 
transverse energy and J/tp production rate on the number iV" par t of partici- 
pant nucleons in high energy heavy ion collisions is extremely important both 
for a better understanding of the initial conditions in the evolution of newly 
created dense matter and becaiise.it. provides the information for discrimi- 
nating among different models BH'EH3#El In this contribution we analyse such 
quantities in the framework of the Dual String Model (DSM). 

We start by building nucleus-nucleus collisions as resulting from superpo- 
sition of nucleon-nucleon collisions, in the way it is done in the Glauber mode] 
approach and its generalisations: in the DSM, i.e., the Dual Parton Modefl 
with the inclusion of strings ,@ the valence quarks of the nucleon produce par- 
ticles, via strings, only once — this is the wounded nucleon model case — and 
production is proportional to the number Na of participant nucleons. As the 
energy and Na increase the role of sea quarks and gluons increases, they in- 
teract and produce, again via strings, particles, and the number of collisions 
v becomes the relevant parameter. One should notice that multiple inelastic 
scattering may occur either internally within a given nucleon-nucleon collision 
or externally involving interactions with different nucleons. 

Following Ref. and taking into account the above basic properties, we 
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now write an expression for the particle pseudo-rapidity density. 
dN 



dy 



= N A [2 + 2{k - l)a] h + {is NA - N A )2kah, (1) 

N A N A 

where h is the height of the valence- valence rapidity plateau, a is the relative 
weight of the sea-sea (including gluons) plateau and k is the average number of 
string pairs per collision. Elementary multi-scattering arguments*! give v Na = 

N^ 3 . However, as we mentioned above, the diagram corresponding to sea-sea 
scattering can be iterated with k > 1 being, in general, a function of energy. 
The number of nucleon-nucleon collisions is, of course, Na + { v n a — ^a) = 
v Na , and the number of strings is N s — Na [2 + 2(k — 1)1 + {v N — N A )2k = 
2ku NA . The first term on the right-hand side of Eq. (|l|) is just a sum over 
nucleon-nucleon scattering contributions (including internal parton multiple 
scattering) and we can thus write 



dN 
dy 



_ dN 

N A N A dy 



Na + {v NA - N A )2kah (2) 

pp 

If external multiple scattering is absent, by putting v Na = Na, one obtains 
the wounded nucleon model limit; if k 3> 1 we obtain the limit in which 
multiple scattering dominates. 

In Fig. y, together with the PHOBOS data, we have presented the quan- 
tity jj- 4?7 as function of the cm. energy y/s for the wounded nucleon 

A y N A N A 

model limit — solid line — and the multiple scattering dominance limit — dotted 
line. Assuming that h and a are energy independent (constant plateaus), the 
energy dependence of dN/dy\ pp , obtained from a parametrisation of exper- 
imental data,Q fixes the energy dependence of k. We find a = 0.05 and 

h = ' 75 - nn 

In the DSM, strings may interact by fusingoEj in the transverse plane 
of interaction thus modifying the number and the distributions of produced 
particles: in particular, due to the vector nature of the colour charge, a cluster 
of m strings will emit fewer particles that m separate strings .Ej 

The number of strings coming from nucleon multiple scattering — the 
second term in Eq. (d) — is Na(N^/ 3 — l)2k and they occupy the transverse 
interaction area S N , which, for central collisions, is approximately given 

by S Na ~ 7r fl.l4/vj/ 3 ^ , such that the dimensionless transverse density 
parameter r) is 



(—) 
U.14/ 



2 2kN l f{N l f-l), (3) 
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Figure 1. Pseudo-rapidity density nor- 
malised per participant pair as a function of 
cm. energy. The lines give predictions for 
the wounded nucleon model (solid line), the 
pure multicollision approach (dotted line), 
and the Dual String Model, without fu- 
sion Eq. (dash-dotted line) and with 
fusion Eq. M) (dashecLJine). AA points 
are tak en from Ref. |s]|b||5|, pp and pp from 

Ref . nq,nji3 13I 



Figure 2. Central charged particle rapidity 
density per participating pair as a function 
of the number of participants. Results of HI- 
JING (histograms), EKRT predictions (dot- 
dashed lines) and DSM predictions (solid 
lines) for central Au+Au collisions at = 
56,130,200 AGeV. Also shown are results 
from pp and pp collisions and PHOBOS data 
(Everything in the figure except the DSM 
curves is taken from Ref. H.) 



where r s ~ 0.2 fm is the string transverse section radius. Note that rj increases 
with Na and, via k, also with \fs. n 
When fusion occurs, Eq. (||) becomesO 



1 dN 



N A N A 



dN 
dy 



F(r))(N\ /3 - l)2kah, 



(4) 



where F(rj) is the particle production reduction factor 



F( V ) 



1 - e- ? > 



V 



(5) 



It can now easily be shownB that the DSM with fusion predicts satura- 
tion of the particle rapidity densities per participant pak_nLnucleons as Na 
increases. This prediction is compared to other modelsEEj£j 
to experimental data from WA980 in figure [| 



in figure y and 
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Figure 3. (a) Charged particle density per participant nucleon versus the number of partic- 
ipants; (b) [absolute charged particle density versus the number of participants. The data 
from WA98EU-efer to 158A GeV Pb+Pb collisions (filled circles), the open circle refers to 
pp collisionsEil; the solid line results from Eq. M) with a = 0.11 and h = 0.77. 



Furthermore, it is to be noted that the predictions for particle densities 
in central Pb+Pb collisions of the DSM without fusion and of the DSM with 
fusion are very different at V* = 200 AGeV (RHIC) and at yfs = 5.5 ATeV 
(LHC) as can be seen in the following table: 



cm. energy 


200 AGeV 


5.5 ATeV 


without fusion 


1500 


4400 


with fusion 


700 


1400 



Of course this model is essentially soft. The parameters of the elementary 
collision densities, h and a, were assumed constant, all the energy dependence 
being attributed to the parameter k, the average number of string pairs per 
elementary collision. If h and a are allowed to grow with energy, as a result, 
for instance, of semi-hard effects, the parameter k may then have a slower 
increase than the one obtained here. 

Finally, one should consider the idea that string fusion eventually leads to 
a situation of percolationoEJ with the formation of extended regions of colour 
freedom, with the features of the expected Quark-Gluon Plasma. Indeed 
the parameter r\ at the CERN-SPS has the value r\ « 1.8, larger than the 
critical density (j] c ~ 1.12 -j- 1.17) which means that percolation transition is 
already taking place at \/s — 20 AGeVp-even allowing for non-uniform matter 
distribution in the nucleus {rj c w 1.5) p3 this result is valid even with k = 1. 
The observed anomalous J/ ip suppressions may then be a signature of the 
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percolation transition to the Quark-Giuon Plasma.E-jli3 

Indeed, in our simple approach,!^ J/ip and Drell-Yan production are 
treated as rare events: this implies that their ratio is given by the product of 
two functions, one describing absorption of J/ip (which we assume as usual to 
be exponential in the amount of matter longitudinally traversed), the other 
describing J/ip (cc) suppression due to Debye screening. If we take the drastic 
position that the latter is 100% effective if there is percolation, and ineffective 
otherwise, then screening is described by the probability of non-percolation, 
which can be parametrised as 

(6) 

with a c a parameter linked to the finite size of the nuclear system. Thus we 
see that the onset of the phase transition is characterised by a change in the 
curvature of the Jji\) over D.Y. ratio from positive (during absorption) to 
negative. This however is only a qualitative description: a quantitative one 
should probably take into account more details of the process (e.g., geometry 
varying with impact parameter, resonances, . . . ). 

In conclusion, the DSM is a model with two components, the valence- 
valence component and the sea-sea component, the sea-sea component increas- 
ing its importance with energy and number of participants. This is somewhat 
similar to the HIJING Monte Carlo model, with soft and hard components. 
On the other hand, with fusion the DSM behaves, for large Na, similarly 
to the EKRT model, but with strict saturation of the particle density per 
participant nucleon. However, in the original EKRT model the saturation 
criterion in the transverse plane is stronger than in case of fusion of strings. 
Here, saturation in the interaction area is asymptotic (when 77 — > 00) while in 
the EKRT model it occurs at finite density. This causes the decrease of the 
particle density with Na in the EKRT original model. 

Probably different explanations, such as the ones based on string fusion, 
parton saturation, parton shadowing, are in some sense dual and refer to 
the same underlying physics .El What is becoming clear is that saturation of 
particle density puts strong constraints in models, and limits the rise of the 
(pseudo-)rapidity plateau at RHIC and LHC. 
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Pn, 



1 + exp 



V~Vc 



References 

1. M. Gyulassy, these proceedings. 

2. B.B. Back et al., PHOBOS Collaboration, Phys. Rev. Lett. 85, 3100 
(2000). 

3. M.M. Aggarwal et al., WA98 Collaboration, preprint nucl-ex/0008004, 
(WA98/CERN). 

4. C. Lourengo, talk given at the Heavy Ion Day meeting in Lisbon, Portu- 
gal, April 2000. 

5. X.-N. Wang and M. Gyulassy, preprint nucl-th/0008014. 

6. J. Dias de Deus and R. Ugoccioni, Phys. Lett. B 491, 253 (2000). 

7. A. Capella, U.P. Sukhatme, C.I. Tan and J. Tran Thanh Van, Physics 
Reports 236, 225 (1994). 

8. N.S. Amelin, M.A. Braun and C. Pajares, Phys. Lett. B 306, 312 (1993); 
N.S. Amelin, M.A. Braun and C. Pajares, Z. Phys. C 63, 507 (1994). 

9. N. Armesto and C. Pajares, Int. J. Mod. Phys. A 15, 2019 (2000). 

10. M. Adamus et al., NA22 Collaboration, Z. Phys. C 37, 215 (1988). 

11. G.J. Alner et al., UA5 Collaboration, Physics Reports 154, 247 (1987). 

12. R.E. Ansorge et al., UA5 Collaboration, Z. Phys. C 43, 357 (1989). 

13. F. Abe et al., CDF Collaboration, Phys. Rev. D 41, 2330 (1990). 

14. N. Armesto, M.A. Braun, E.G. Ferreiro and C. Pajares, Phys. Rev. Lett. 
77, 3736 (1996); M. Nardi and H. Satz, Phys. Lett. B 442, 14 (1998). 

15. J. Dias de Deus, R. Ugoccioni and A. Rodrigues, Eur. Phys. J. C 16, 537 
(2000). 

16. T.S. Biro, H.B. Nielsen and J. Knoll, Nucl. Phys. B 245, 449 (1984). 

17. M.A. Braun and C. Pajares, Eur. Phys. J. C 16, 349 (2000). 

18. J. Dias de Deus and R. Ugoccioni, Phys. Lett. B 494, 53 (2000). 

19. K.J. Eskola, K. Kajantie, P.V. Ruuskanen and K. Tuominen, Nucl. Phys. 
B 570, 379 (2000). 

20. K.J. Eskola, K. Kajantie and K. Tuominen, preprint JYFL-3700 and HIP- 
2000-45/TH (hep-ph/0009246), Univ. of Jyvaskyla and Univ. of Helsinki. 

21. C. De Marzo et al., Phys. Rev. D 26, 1019 (1982). 

22. A. Rodrigues, R. Ugoccioni and J. Dias de Deus, Phys. Lett. B 458, 402 
(1999). 

23. M.C. Abreu et al., NA50 Collaboration, Phys. Lett. B 477, 28 (2000). 



ugoccioni: submitted to World Scientific on February 1, 2008 6 



